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First measurements of charge-independent correlations on angular difference variables r\\ — V2 
(pseudorapidity) and <j>\ — <p2 (azimuth) are presented for primary charged hadrons with transverse 
momentum 0.15 < p t < 2 GeV/c and \r]\ < 1.3 from Au-Au collisions at ^/snn = 130 GeV. Strong 
charge-independent angular correlations are observed associated with jet-like structures and elliptic 
flow. The width of the jet-like peak on r/i — r/2 increases by a factor 2.3 from peripheral to central 
collisions, suggesting strong coupling of semi-hard scattered partons to a longitudinally-expanding 
medium. New methods of jet analysis introduced here provide evidence for nonperturbative QCD 
medium effects in heavy ion collisions. 
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PACS numbers: 24.60.-k, 24.60.Ky, 25.75.Gz 



I. INTRODUCTION 

Correlations and fluctuations can provide essential in- 
formation on the nature of the medium produced in ultra- 
relativistic heavy ion collisions |iJ. 121 L| In-medium mod- 
ification of parton scattering and fragmentation of ener- 
getic partons and the bulk medium can be observed via 
large-momentum-scale correlations. Charge-independent 
angular correlations result from initial-state multiple 
scattering (Cronin effect @|, hard parton scattering yjj]) 
with subsequent in-medium parton dissipation and el- 
liptic flow. Medium modification of minimum-bias, semi- 
hard parton scattering and fragmentation (i.e., minijets) 
is the subject of this paper. 

Previous studies of parton-medium interactions have 
included angular correlations of high-p t particles based 
on a leading-particle technique (e.g., leading-particle pt > 
4 GeV/c, associated particle pt < 4 GeV/c) in which the 
away-side jet structure was observed to be strongly re- 
duced in central Au-Au collisions 0. Theoretical de- 
scriptions of parton energy loss and medium-modified 
fragmentation include perturbative quantum chromody- 
namics (pQCD) based jet-quenching models 0, and 
parton recombination models |9(. 

In this Letter we report the first measurements in 
heavy ion collisions of charge-independent joint autocor- 
relations [lOj on angular difference variables 4>a = 4>\~<\>2 
(azimuth) and t]a = f?i — 772 (pseudorapidity) for charged 
particles with 0.15 < pt < 2 GeV/c. This analysis is 
based on s/s^n = 130 GeV Au-Au collisions observed 
with the STAR detector ^1 at the Relativistic Heavy 
Ion Collider (RHIC). The measurements involve no a 
priori jet model (specifically, no trigger particle). The 
autocorrelation technique, combined with the large an- 
gular acceptance of the STAR detector, enables statisti- 
cally weak correlation structures, which are individually 
undetectable but occur multiple times in each event, to 
be measured in the aggregate with good statistical ac- 
curacy. These novel low-p t measurements reveal jet-like 
correlations which suggest that in central Au-Au colli- 
sions strong non-perturbative coupling of partons to a 
longitudinally-expanding medium jl2j produces dramatic 
changes in the angular distribution of parton fragments 
(hadrons) not anticipated by pQCD-based theory [EL Il3|. 



considered in a separate analysis |14|. We further project 
the 4D (four dimensional) subspace onto a 2D subspace 
of angular difference variables to form a joint autocorre- 
lation which is observed to retain essentially all angular 
correlation information The autocorrelations ob- 

tained in this study thus access the complete CI angular 
structure of two-particle density p(p~l,p~2). 

Differential correlation analysis is achieved by compar- 
ing an object distribution p S ib of particle pairs taken from 
single events (sibling pairs) with a reference distribution 
Pmix where each particle in the pair is taken from differ- 
ent but similar events (mixed pairs). The corresponding 
correlation function and pair-number density ratio are 
defined by 

C(p\ , P2 ) = Psib (fil, P2) - Pmix (Pi , P2 ) 
r(Pl,P2) = Pstb(Pl,P2)/ Prmx(Pl,P2)- (1) 

Pair densities p(pi,P2) are first projected onto subspaces 
(771,772), (01,02) as histograms n ab ~ e x e y p(x a , y b ), 
where ab are 2D bin indices and e x , e y are bin widths 
on x,y S {77, <f>, 77A, 4>a}- Sibling- and mixed-pair his- 
tograms are separately normalized to the total number 
of detected pairs in each event class: h a b — n a b/ J2 a b Uab - 
Normalized pair- number ratios f a b — n a b yS ib/n a b,mix are 
the basis for this analysis. Ratios are formed from sub- 
sets of events with similar centrality (multiplicities differ 
by < 100, except < 50 for most-central) and primary- 
vertex location (within 7.5 cm along the beam axis) and 
combined as weighted (by sibling pair number) averages 
within each centrality class. If correlation structure is in- 
variant on sum variables 771+772 and 0i+02 (stationarity), 
as observed previously in heavy ion collisions [15j, his- 
tograms f a b may be projected by averaging along parallel 
diagonals to form 2D joint autocorrelations fwi on 
difference variables 77a,</>a- Autocorrelation details are 
described in [lfl. Il7| . The autocorrelations in this study 
should not be confused with conditional distributions ob- 
tained from leading-particle jet analyses. 2D joint auto- 
correlations compactly represent all angular correlations 
on 4D subspace (771, 772, 4>i, 4>2) without information loss 
(distortion). 

III. DATA 



II. ANALYSIS METHOD 

Our goal is to access the complete charge-independent 
(CI - all charged particles) structure of two-particle den- 
sity p(pitP2.) without imposing a correlation model. The 
full two-particle momentum space is projected onto an- 
gular subspace (771, 772, <fii, </>2), integrating over a specific 
transverse momentum interval. Correlation structure on 
transverse momentum with specific angular constraints is 



Data for this analysis were obtained with the STAR 
detector ^l| using a 0.25 T uniform magnetic field par- 
allel to the beam axis. Event triggering and charged- 
particle measurements with the time projection chamber 
(TPC) are described in Track definitions, track- 

ing efficiencies, quality cuts and primary-particle defini- 
tion are described in 0, 0] . Tracks were accepted in 
\rj\ < 1.3, 0.15 < p t < 2 GeV/c and full azimuth. Parti- 
cle identification was not implemented. Corrections were 
made to f for two-track inefficiencies due to track merging 



(a) (b) (a) (b) 




FIG. 1: Perspective views of two-particle CI joint autocorre- 
lations N(r — 1) on (?7a,0a) for central (a) to peripheral (d) 
collisions. 



and intersectin g tr ajectories reconstructed as > 2 parti- 
cles (splitting) |2£j. Small-scale momentum correlations 
due to HBT and Coulomb effects [21| were suppressed by 
eliminating sibling and mixed track pairs with |?7a| < 0.3, 
\M < tt/6, \pti-pa\ < 0.15 GeV/c, if p t < 0.8 GeV/c 
for either particle. These pair cuts have negligible effect 
on the correlations studied here. Four centrality classes 
for 300k events labeled (a) - (d) for central to peripheral 
were defined by cuts on TPC track multiplicity N within 
the acceptance relative to minimum-bias event multiplic- 
ity frequency distribution end-point A^o |22| . which cor- 
responds to the maximum participant number |l8l l23| . 
The four centrality classes used here were defined by 
(d) 0.03 < N/N < 0.21, (c) 0.21 < N/N < 0.56, (b) 
0.56 < N/Nq < 0.79 and (a) N/N > 0.79, correspond- 
ing respectively to fraction of total cross section ranges 
40-70%, 17-40%, 5-17% and 0-5%. 



IV. TWO-PARTICLE DISTRIBUTIONS 

Plotted in Fig.^are perspective views of CI joint auto- 
correlations of quantity N(f — 1) (measuring the density 
of correlated pairs per final-state particle 15], typically 
O(l) for all centralities, N is the mean multiplicity in the 
acceptance) for four centrality classes of Au-Au collisions. 
N(f— 1) would be independent of centrality if Au-Au col- 
lisions were linear superpositions of p-p collisions (partic- 
ipant scaling). The distributions in Fig.^are dominated 
by 1) a ID quadrupole component oc cos(20a) conven- 
tionally attributed to elliptic flow; 2) a ID dipole com- 



FIG. 2: The same data as in Fig.0 but with 7)A-independent 
dipole and quadrupole components on 0a (see text) sub- 
tracted to reveal 'same-side' (\4>a\ < tt /2) structures which 
can be associated with minijets. 

ponent oc cos(</)a) associated with transverse momentum 
conservation in a thermal system, and 3) a 2D 'same-side' 
(|<M| < 7r /2) peak. The same-side peak is assumed to be 
associated with parton fragmentation to hadrons, albeit 
for fragments with much lower p t s than are considered in 
a conventional jet analysis. Since no leading or trigger 
particle is invoked in this model-independent analysis we 
refer to the corresponding fragment systems as minijets 
(minimum-bias jets). 

Momentum conservation in a thermal multiparticle 
system should result in a dipole cos(0a) angular cor- 
relation [24|]. We also expect back- to-back or away-side 
(4>a ~ tt) azimuth correlations from momentum conser- 
vation in parton scattering (dijets). However, at low p t 
the away-side jet structure is broad, and indistinguish- 
able from the dipole cos{4>&) component describing mo- 
mentum conservation in the bulk system. We subtract 
dipole and quadrupole cos (2(/>a) components from dis- 
tributions in Fig. ^ to obtain Fig. [5] by minimizing 77A- 
independent sinusoidal residuals on the away side region 
(\<Pa\ > tt/2) and for |?7a| ~ 2. The small excess in the 
(0,0) bins is due to conversion-electron pair contamina- 
tion. The same-side 2D peaks in this figure are the main 
subject of this analysis. 

We observe that the away-side region in Fig- El is fea- 
tureless, even for the most peripheral collisions. If Lund- 
model strings |25| remained dynamically relevant in the 
final stage of heavy ion collisions we would expect, in 
the accepted pt interval, significant correlation structure 
on the away side of Fig. [21 a prominent ID gaussian on 



5 



t]a approximately symmetric on azimuth and due to lo- 
cal charge conservation on z (during longitudinal string 
fragmentation) and coupling of z to 77 due to longitudinal 
Hubble expansion, as observed in p-p collisions 
The absence of such structure suggests that longitudinal 
strings play no significant role in the final stage of Au-Au 
collisions, even for the most peripheral collisions in this 
study. That trend is consistent with the centrality depen- 
dence of net-charge correlations in which structure char- 
acteristic of string fragmentation is strongly supp ressed 
with increasing centrality of Au-Au collisions |l5| . 

The same-side peak isolated in Au-Au collisions by 
the multipole subtraction varies strongly with central- 
ity, transitioning from significant elongation on azimuth 
difference a for p-p collisions (2^ to dramatic broaden- 
ing along ?7a for the more central Au-Au collisions (note 
the non-unit aspect ratio of these 2D plots). HBT and 
Coulomb pair cuts (see Sec. 1111(1 reduce the bins near- 
est (0,0) by 10% or less. ID projections of data and 2D 
model fits (discussed below) onto difference variables a 
and 77A are shown in Fig. [3] Solid dots and curves (open 
triangles and dashed curves) correspond to 77 a (0a) pro- 
jections. 



V. ERRORS 

Statistical errors for joint autocorrelations approxi- 
mately double as |tja| increases from to 2 because of 
the bounded rj acceptance, but are uniform on 0a be- 
cause the azimuthal acceptance is continuous (periodic) 
in STAR. Statistical errors for f at |?7a| =0 vary from 
0.0001 for central collisions to 0.001 for peripheral colli- 
sions. Statistical errors for N(f — 1) (~ 0.1) are nearly 
independent of centrality. Systematic errors were esti- 
mated as in |l8j . Contamination from photon conversions 
to pairs is significant only within the bin defined by 
1 77A I < 0.1, |0a| < 0.1 which was excluded from model 
fits. The dominant source of systematic error is non- 
primary background, mainly weak-decay daughters |19| . 
whose correlation with primary particles is unknown and 
is estimated by assuming those correlations vary from 
zero to the measured correlation amplitude for primary 
particles • Total systematic errors for data presented 
in Fig. n are ±7% of signal, but increase to ±8% for 
|?7a| < 0.5 and to ±11% for |0 A | < 0.05. Correlations 
from resonance (p°,oj) decays are ~3% of the peaks at 
(0,0) in Fig.Hin |r7 A | < 0.5, |0 A | < 2 H|. 



VI. MODEL FITS 

Joint autocorrelations, as in Fig. ^ but without factor 
JV, were fitted with a model function consisting of dipole 
and quadrupole terms on 0a, a ID gaussian on 77A and a 
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FIG. 3: Projections of 2D CI joint autocorrelations N(r — 1) 
in Fig.|2|onto difference variables 77a (solid dots) and 0a (open 
triangles) . The solid (dashed) curves represent corresponding 
projections of 2D analytical model fits to the data. 2D peaks 
are substantially reduced in amplitude by ID projections. 



2D same-side gaussian on (77A, 0a), plus a constant offset 
F = A 0A cos(0 A ) + A 2 < l>& cos(2 A ) + A e \^°o > 

(2) 
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Best-fit parameters for the model fits shown in Fig.[3]arc 
listed in Table [fl including mean multiplicity factor N 
as in Figs. Ill il plus tracking efficiency correction factor 
S |2^. Those fit parameters confirm that with increas- 
ing centrality 2D peak structures exhibit 1) strong and 
non-monotonic amplitude variation, 2) strong 77A width 
increase and 3) significant 0a width reduction. 



TABLE I: Parameters and fitting errors (only) for model fits 
[Eq. (J2J] to joint autocorrelation data in Fig. Q for centrality 
bins (a) - (d) (central - peripheral). Total systematic error 
for efficiency-corrected amplitudes is 11% |2flj . 



centrality 


(d) 


(c) 


(b) 


(a) 


error" (%) 


5[29] 

N 


1.19 
115.5 


1.22 
424.9 


1.25 
790.2 


1.27 
983.0 


8 (syst) 


SNAi 


1.93 
0.61 
0.58 


3.23 
0.55 
1.05 


3.72 
0.54 
1.34 


3.10 
0.53 
1.36 


5-2 

4- 2 

5- 2 


SNA 
cro 


0.60 
1.11 


0.32 
0.24 






0.16-0.1 6 
28-22 


SNA 2 


-0.67 


-0.55 


-0.67 


-0.58 


C 


SNA^ & 
SNA 24>& 


-0.31 
1.05 


-0.76 
2.72 


-0.97 
1.30 


-0.74 
0.32 


22-5 
2-17 


X 2 /DcF 


439 
31 <; 


419 
31fi 


675 
31(5 


415 
31 (i 





"Range of fitting errors in percent from peripheral to central. 
''Magnitude of fitting errors. 
c Fixed by normalization of f 
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VII. DISCUSSION 

In Fig. 01 the same-side peak amplitudes and widths 
from model fits are plotted vs centrality measure v |23| 
(mean participant path length in terms of the number of 
encountered nucleons, ~ 2Nbi n /N part ) along with mea- 
surements obtained from p-p collision data 27}. The 
same-side peaks in Fig. differ strongly from those 
for p-p collisions, where for the latter a 2D gaussian 
peak significantly elongated on azimuth dominates the 
same-side structure, widths on 77 a and (j>A being ~ 0.5 
and 0.7 respectively |23|- The similar same-side peak 
widths for mid-peripheral Au-Au collisions in this anal- 
ysis (panel (d), v ~ 2.5) are consistent with the p-p re- 
sult. In central Au-Au collisions however, the widths of 
the same-side peak reverse the sense of the asymmetry: 
the peak is dramatically elongated on ?7a, the width ratio 
increasing to 2.6. 




v 



FIG. 4: Left: Fitted amplitudes and volumes for peaks in 
Fig. |3 plotted vs mean participant path length v from 
TableQ Right: Fitted widths a VA (dots) and (J^ A (triangles). 
Hatched regions show p-p values. Curves guide the eye. 

Same-side peak amplitudes SNA\ (measuring correla- 
tions per final-state particle) in the left panel of Fig. 0] 
increase nearly linearly with path-length as expected for 
independent binary collisions. However, peak volume 
= SNAi a nA <70 A (oc minijet fragment number) has a 
more complex variation, strongly departing from linear v 
scaling (dotted line) above v = 2.5. The dashed curve in 
the left panel is derived from the curves describing am- 
plitude and peak widths. The volume excess beyond the 
linear extrapolation may indicate the onset of a strongly 
dissipative medium in which more correlated fragments 
with less pt result from each scattered parton. The lat- 
ter increase is very likely a lower-p t manifestation of the 
observed suppression of the high-popart of the p t spec- 
trum measured by quantity Raa |3_Qj . It is notable that 
the peak amplitude does not deviate from a linear trend, 
except for the most central point. 

We speculate that the mechanism modifying the same- 
side peak in central Au-Au collisions is strong coupling 
of minimum-bias energetic partons (no high-p t trigger is 
imposed) to a longitudinally- expanding colored medium 
developed in the more central Au-Au collisions. Hadrons 
from fragmenting partons sample the local velocity struc- 
ture of the pre-hadronic parton-medium coupled system. 



Growth of the colored medium and its coupling to the 
fragmenting parton with increasing collision centrality is 
then indicated by increased width on r/& of the same-side 
correlation peak. 

The perturbative QCD expectation for angular corre- 
lations about the jet thrust axis in p-p collisions corre- 
sponds to a nearly-symmetric near-side peak on 4>a). 
That 'in vacuo' result is indeed observed in p-p collisions 
for higher-p t fragments (> 2.5 GeV/c). However, in a 
minimum-bias autocorrelation analysis of p-p data |27j 
strong deviations from expected pQCD angular symme- 
try about the jet thrust axis are observed. The Hijing 
Monte Carlo collision model [f| includes a conventional 
pQCD model of jet production and quenching in A-A 
collisions. The default Hijing same-side peak is observed 
to be symmetric, and the widths on tja and (j>A both in- 
crease by only 10% when jet quenching is imposed [3l| . 
seriously underpredicting the large width increase on r/A 
and contradicting the width decrease on 0a observed in 
the present analysis of Au-Au data. The pQCD jet- 
quenching mechanism in Hijing cannot produce an asym- 
metry on (r/A, <M)j given the symmetry about the parton 
momentum of its perturbative bremsstrahlung quench- 
ing mechanism. Prominent \ow-p t longitudinal string- 
fragment correlations on r/A are observed for all Hijing 
centralities, also contradicting results of the present anal- 
ysis noted above. RQMD [jjj CI correlations are feature- 
less except for small flow- related correlations on 4>a- 

Recently, effects of a flowing medium on parton en- 
ergy loss and fragmentation have been explored theo- 
retically |8|. The premise of that study is that gluon 
bremsstrahlung from energetic partons transiting a col- 
ored medium should be sensitive to the local structure 
of the velocity field on the medium. The model consid- 
ered is uniform medium flow ('directed flow') transverse 
to the energetic parton momentum. A static medium is 
expected to broaden the bremsstrahlung angular distri- 
bution and hence the near-side peak (symmetric broad- 
ening from jet quenching is observed already with Hi- 
jing [3l|]). Medium flow transverse to the parton direc- 
tion was found in @ to shift and distort the fragment- 
energy angular distribution relative to the thrust axis. In 
the LHC context, for 100 GeV jets with typical energy 
angular width ~ 0.05, the effect of the flowing medium 
on the angular distribution was found to be comparable 
to the width magnitudes ('marked medium-induced de- 
viations'). However, the absolute angular changes were 
small. 

In the RHIC context a comparison was made with a 
STAR leading-particle analysis of jet correlations [32| . 
The prediction of § for trigger particles with p t £ [4,6] 
GeV/c is width variation from peripheral to central 200 
GeV Au-Au collisions of 0.35 (symmetric) to 0.4 on az- 
imuth and to 0.56 on pseudorapidity. Those width in- 
creases are similar in magnitude to the symmetric Hijing 
width increases noted above. However, they differ qual- 
itatively from the width decrease from 0.7 to 0.5 on az- 
imuth and the dramatic width increase from 0.5 to 1.4 
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on pseudorapidity observed in the present minimum-bias 
jet study. The study of directed (vector) flow and par- 
ton bremsstrahlung in || does not address the issue of 
longitudinal Hubble (Bjorken) expansion, a tensor aspect 
of the velocity field. Coupling of parton fragmentation 
to the velocity field may be much stronger than what 
can be modelled perturbatively, requiring a nonperturba- 
tive treatment. The analysis in [8| also does not address 
the centrality dependence of angular deformation, which 
is strongly nonlinear on path length as demonstrated in 
Fig. 0] (right panel) . 

Fluctuation analysis has been generally advocated as a 
probe of heavy ion collisions. It is important to note that 
charge-independent number fluctuations observed within 
a given detector acceptance integrate over that accep- 
tance the CI joint autocorrelations presented in this pa- 
per (within a constant offset), as described in The 
constant offset in the joint autocorrelations, dominated 
by participant ('volume') fluctuations, is easily separated 
from the differential angular structure which reveals de- 
tails of nuclear collision dynamics. 

VIII. SUMMARY 

In conclusion, we have for the first time measured 
charge-independent joint autocorrelations on difference 
variables </>a and t/a for Au+Au collisions at ^/sjvat 
= 130 GeV. Low-p t string-fragment correlations which 
appear prominently in p-p collisions are not observed 
for any centrality in this study: longitudinal string de- 
grees of freedom are apparently strongly suppressed even 
for fairly peripheral Au-Au collisions. Other correlation 



structures are observed to have substantial amplitudes. 
In addition to azimuth structures associated with ellip- 
tic flow and transverse momentum conservation we ob- 
serve a near-side peak structure varying from a nearly- 
symmetric shape on (^Ai^a) in peripheral collisions to 
a shape strongly elongated on 77 a in central collisions. 
We interpret the same-side peak as resulting from frag- 
mentation of minimum-bias partons observed with no 
trigger condition (minijets). The trend of minijet an- 
gular deformation, observed in this first jet analysis with 
low-pt hadrons, can be interpreted as a transition from 
in vacuo jet fragmentation in p-p and peripheral Au-Au 
collisions to strong coupling of minimum-bias partons to 
a longitudinally-expanding colored medium in the more 
central collisions as part of a parton dissipation process. 
The concept of parton energy loss in heavy ion colli- 
sions is thereby extended to strongly nonperturbative as- 
pects. Detailed comparisons between data and pQCD- 
based theory reveal several qualitative differences. 

We thank the RHIC Operations Group and RCF at 
BNL, and the NERSC Center at LBNL for their sup- 
port. This work was supported in part by the HENP 
Divisions of the Office of Science of the U.S. DOE; the 
U.S. NSF; the BMBF of Germany; IN2P3, RA, RPL, 
and EMN of France; EPSRC of the United Kingdom; 
FAPESP of Brazil; the Russian Ministry of Science and 
Technology; the Ministry of Education and the NNSFC 
of China; IRP and GA of the Czech Republic, FOM of 
the Netherlands, DAE, DST, and CSIR of the Govern- 
ment of India; Swiss NSF; the Polish State Committee 
for Scientific Research; STAA of Slovakia, and the Korea 
Sci. & Eng. Foundation. 



[1] R. Stock, Nucl. Phys. A661, 282c (1999); H. Heiselberg, 

Phys. Rep. 351, 161 (2001). 
[2] A. Dumitru, R. Pisarski, Phys. Lett. B504, 282 (2001). 
[3] L. M. Bettencourt, K. Rajagopal and J. V. Steele, Nucl. 

Phys. A693, 825 (2001). 
[4] M. Gazdzicki, A. Leonidov, G. Roland, Eur. Phys. J. C6, 

365 (1999). 

[5] X.-N. Wang, M. Gyulassy, Phys. Rev. D 44, 3501 (1991). 
[6] Q. Liu and T. A. Trainor, Phys. Lett. B567, 184 (2003). 
[7] C. Adler et al, Phys. Rev. Lett. 90, 082302 (2003). 
[8] N. Armesto, C. A. Salgado and U. A. Wiedemann, Phys. 

Rev. Lett. 93, 242301 (2004) 
[9] R.Hwa and C.Yang, Phys. Rev. C 66, 025205 (2002). 
[10] An autocorrelation is a projection by averaging from sub- 
space (x\,X2) onto difference variable xa = x\ — X2- A 
joint autocorrelation is a simultaneous projection onto 
two difference variables. 
[11] K. H. Ackermann et al, Nucl. Instrum. Meth. A 499, 

624 (2003); see other STAR papers in volume A499. 
[12] H. R. Schmidt and J. Schukraft, J. Phys. G: Nucl. Part. 

Phys. 19, 1705 (1993). 
[13] H. Sorge, H. Stocker, W. Greiner, Nucl. Phys. A498, 
567c (1989); Ann. Phys. (N.Y.) 192, 266 (1989). 



[14] J. Adams et al. (STAR Collaboration) , nucl-ex/0408012 
[15] J. Adams et al. (STAR Collaboration), nucl-ex/0406035 
[16] T. A. Trainor, R. J. Porter and D. J. Prindle, J. Phys. G 
31, 809 (2005). 

[17] D. J. Prindle and T. A. Trainor, hep-ph/0506173 To ap- 
pear in: Proceedings of the MIT Workshop on Correla- 
tions and Fluctuations in Relativistic Nuclear Collisions, 
Cambridge, Massachusetts, 21-23 April, 2005. 

[18] J. Adams et al. (STAR Collaboration), Phys. Rev. C 71 
064906 (2005). 

[19] C. Adler et al, Phys. Rev. Lett. 87, 112303 (2001); ibid. 
89, 202301 (2002). 

[20] Track cuts required minimum track separations 5, 12 and 
20 cm at radii 50, 127 and 200 cm from the TPC axis. 
Crossing pairs with separations less than 10 cm (z) and 
30 cm (azimuth) at mid-radius were also excluded. 

[21] C. Adler et al. Phys. Rev. Lett. 87, 082301 (2001). 

[22] No, the half-max point at the end of the minimum-bias 
distribution plotted as do / dN^[ A , is an estimator on mul- 
tiplicity N for the maximum number of participant nu- 
cleons; N/N ~ N part /Np a rt,max within 4%. 

[23] v estimates the mean participant path length as a 
number of encountered nucleons. For this analysis 



8 



v = 5.5(iV/iVo) 1/3 ^ 5.5(N part /N part ,m ax ) 1/:i ~ 

2Nbi n /Npart, based on Glauber-model simulations. N par t 

is the number of participants and Nu n is the number of 

binary collisions. 
[24] H. Feshbach, A. Gal and J. Hiifner, Ann. Phys. (N.Y.) 

66, 20 (1971) (see Sec. 3.2). 
[25] B. Andersson, G. Gustafson, G. Ingelman and T. 

Sjostrand, Phys. Rep. 97, 31 (1983). 
[26] J. Whitmore, Phys. Rep. 27, 187 (1976). 
[271 R- J- Porter and T. A. Trainor (STAR Collaboration), 

hep-ph/0406330 



[28] R. Ray and R. Longacre, nucl-ex/0008009 

[29] Extrapolation factors S for N(f — 1) |18l | correct for con- 
tamination and inefficiency 19]. Systematic error in S 
was estimated to be ±8%. Total systematic error for ex- 
trapolated quantities in Table[I]was 11%. 

[30] J. Adams et al, Phys. Rev. Lett. 91, 172302 (2003). 

[31] Q. J. Liu, D. J. Prindle and T. A. Trainor, to be published 
in Phys. Lett B (2005), hep-ph/0410180 ). 

[32] F. Wang (STAR Collaboration), J. Phys. G 30, S1299 
(2004). 



